Draft version December 2, 2010 

Preprint typeset using I^T^jK style emulateapj v. 08/29/06 



GALAXY CLUSTERS SELECTED WITH THE SUNYAEV-ZEL'DOVICH EFFECT 
FROM 2008 SOUTH POLE TELESCOPE OBSERVATIONS 

K. Vanderlinde,^ T. M. Crawford, T. de Haan/ J. P. Dudley/ L. Siiaw/ " P. A. R. Ade,"^ K. A. Aird,'' 
B. A. Benson/'^'** L. E. Bleem.^" M. Brodwin/" J. E. Carlstrom,^-'*'*-^ C. L. Chang, ^''^ A. T. Crites,^'"' S. Desai," 
M. A. DoBBS,^ R. J. FoLEY,^" E. M. George,^ M. D. Gladders,^'-' N. R. Hall,^^ w. Halverson,^'' F. W. High,!" 
G. P. Holder,^ W. L. Holzapfel,' J. D. Hrubes,^ M. Joy,^-' R. Keisler,^-^ L. Knox," A. T. Lee,"^ " E. M. Leitch,^'^ 

A. Loeiir,"' M. Lueker,^ D. P. Marrone,^'*' J. J. McMahon,^ J. Meiil,^-'' S. S. Meyer,^'**--' -' J. J. Moiir,^''^"-^" 
T. E. Montroy,^^ C.-C. Ngeow,"'^^ S. Padin,^-^ T. Plagge,^-^ C. Pryke,^'^-** C. L. Reiciiardt,^ A. Rest," J. Ruel," 
J. E. Ruhl,^^ K. K. Schaffer,^** E. Siiirokoff,^ J. Song," H. G. Spieler, B. Stalder,"' Z. Staniszewski,^^ 
A. A. Stark, ^° C. W. Stubbs,"'^" A. van Engelen,^ J. D. Vieira,^''' R. Williamson,^'^ Y. Yang," O. Zahn,'' and 

A. Zenteno^^'^^ 
Draft version December 2, 2010 

ABSTRACT 

We present a detection-significance-limited catalog of 21 Sunyaev-Zel'dovich-selected galaxy clusters. 
These clusters, along with 1 unconfirmed candidate, were identified in 178 deg^ of sky surveyed in 
2008 by the South Pole Telescope to a depth of 18 /iK-arcmin at 150 GHz. Optical imaging from the 
Blanco Cosmology Survey (BCS) and Magellan telescopes provided photometric (and in some cases 
spectroscopic) redshift estimates, with catalog redshifts ranging from z = 0.15 to z > 1, with a median 
z = 0.74. Of the 21 confirmed galaxy clusters, three were previously identified as Abell clusters, three 
were presented as SPT discoveries in Staniszewski et al. (2009), and three were first identified in a 
recent analysis of BCS data by Menantcau ct al. (2010); the remaining 12 clusters arc presented for 
the first time in this work. Simulated observations of the SPT fields predict the sample to be nearly 
100% complete above a mass threshold of M200 ~ 5 x IO^^Mq/i"^ at z = 0.6. This completeness 
threshold pushes to lower mass with increasing redshift, dropping to 4 x IQ^^ MQh~^ at z = 1. 
The size and redshift distribution of this catalog are in good agreement with expectations based on 
our current understanding of galaxy clusters and cosmology. In combination with other cosmological 
probes, we use this cluster catalog to improve estimates of cosmological parameters. Assuming a 
standard spatially flat wCDM cosmological model, the addition of our catalog to the WMAP 7-year 
results yields as = 0.81 ± 0.09 and w = — 1.07 ± 0.29, a ~ 50% improvement in precision on both 
parameters over WMAP7 alone. 

Subject headings: galaxies: clusters: individual, cosmology: observations 
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1. INTRODUCTION 

The most massive dark matter halos to have formed 
so far have characteristic masses of 10^^ to 10^'' solar 
masses. Although dark matter makes up the vast major- 
ity of the mass of these objects, most observation signa- 
tures result from baryons. A small fraction of the baryons 
in these massive halos eventually cools to form stars and 
galaxies, and it was through the light from these galaxies 
that the most massive halos were first identified. Because 
of this, we generally refer to these objects as galaxy clus- 
ters, despite the small contribution of galaxies to their 
total mass. 

Galaxy clusters are tracers of the highest peaks in the 
matter density field and, as such, their abundance is ex- 
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ponentially sensitive to the growth of structure over cos- 
mic time. A measurement of the abundance of galaxy 
clusters as a function of mass and redshift has the power 
to constrain cosmological parameters to unprecedented 
levels (Wang & Steinhardt 1998; Haiman et al. 2001; 
Holder et al. 2001; Battye & Weller 2003; Molnar et al. 
2004; Wang ct al. 2004; Lima & Hu 2007), assuming 
that the selection criteria are well understood. To use- 
fully constrain the growth history of large scale struc- 
ture, a sample of galaxy clusters must cover a wide red- 
shift range. Furthermore, the observable property with 
which the clusters are selected should correlate strongly 
with halo mass, which is the fundamental quantity pre- 
dicted from theory and simulations. The thermal Sun- 
yaev Zel'dovich (SZ; Sunyaev & Zel'dovich 1972) signa- 
tures of galaxy clusters provide nearly these selection 
criteria. Surveys of galaxy clusters based on the SZ ef- 
fect have consequently been eagerly anticipated for over 
a decade. This paper presents the first cosmologically 
meaningful catalog of galaxy clusters selected via the 
thermal SZ effect. 

1.1. The Thermal SZ Effect 

The vast majority of known galaxy clusters have been 
identified by their optical properties or from their X-ray 
emission. Clusters of galaxies contain anywhere from 
several tens to many hundreds of galaxies, but these 
galaxies account for a small fraction of the total bary- 
onic mass in a cluster. Most of the baryons in clus- 
ters are contained in the intra-cluster medium (ICM), 
the hot 10^ — 10^ K) X-ray-emitting plasma that per- 
vades cluster environments. Sunyaev & Zel'dovich (1972) 
noted that this same plasma should also interact with 
cosmic microwave background (CMB) photons via in- 
verse Compton scattering, causing a small spectral dis- 
tortion of the CMB along the line of sight to a cluster. 

The thermal SZ effect has been observed in dozens of 
known clusters (clusters previously identified in the opti- 
cal or X-ray) over the last few decades (Birkinshaw 1999; 
Carlstrom et al. 2002). However, it was not until very 
recently that the first previously unknown clusters were 
identified through their thermal SZ effect (Staniszewski 
et al. 2009). This is mostly due to the small amplitude 
of the effect. The magnitude of the temperature distor- 
tion at a given position on the sky is proportional to the 
integrated electron pressure along the line of sight. At 
the position of a massive galaxy cluster, this fluctuation 
is only on the order of a part in 10**, or a few hundred 
liK^. It is only with the current generation of large (~ 
kilopixel) detector arrays on 6-12 m telescopes (Fowler 
et al. 2007; Carlstrom et al. 2009) that large areas of sky 
are being surveyed to depths sufficient to detect signals 
of this amplitude. 

A key feature of the SZ effect is that the SZ surface 
brightness is insensitive to the redshift of the cluster. As 
a spectral distortion of the CMB (rather than an intrin- 
sic emission feature), SZ signals redshift along with the 
CMB. A given parcel of gas will imprint the same spec- 

^ Throughout this work, the unit K refers to equivalent flue- 
tuations in the CMB temperature, i.e., the level of temperature 
fluctuation of a 2.73 K blackbody that would be required to pro- 
duce the same power fluctuation. The conversion factor is given by 
the derivative of the blackbody spectrum, 4^, evaluated at 2.73 K. 



tral distortion on the CMB regardless of its cosmologi- 
cal redshift, depending only on the electron density 
and temperature Tg. This makes SZ surveys an excellent 
technique for discovering clusters over a wide redshift 
range. 

Another aspect of the thermal SZ effect that makes it 
especially attractive for cluster surveys is that the inte- 
grated thermal SZ flux is a direct measure of the total 
thermal energy of the ICM. The SZ flux is thus expected 
to be a robust proxy for total cluster mass (Barbosa et al. 
1996; Holder & Carlstrom 2001; Motl et al. 2005). 

A mass-limited cluster survey across a wide redshift 
range provides a growth-based test of dark energy to 
complement the distance based tests provided by su- 
pernovae (Perlmutter et al. 1999; Schmidt et al. 1998). 
Recent results (e.g., Vikhlinin et al. 2009; Mantz et al. 
2010), have demonstrated the power of such tests to con- 
strain cosmological models and parameters. 

1.2. The SPT SZ Cluster Survey 

The South Pole Telescope (SPT) (Carlstrom et al. 
2009) is a 10-meter off-axis telescope optimized for 
arcmimite-rcsolution studies of the microwave sky. It is 
currently conducting a survey of a large fraction of the 
southern sky with the principal aim of detecting galaxy 
clusters via the SZ effect. In 2008, the SPT surveyed 
~ 200 deg^ of the microwave sky with an array of 960 
bolometers operating at 95, 150, and 220 GHz. Using 
40 deg^of these data (and a small amount of overlap- 
ping data from 2007), Staniszewski et al. (2009) (here- 
after S09) presented the first discovery of previously un- 
known clusters by their SZ signature. Lueker et al. (2010) 
(hereafter LIO) used ~ 100 deg^ of the 2008 survey 
to measure the power spectrum of small scale temper- 
ature anisotropies in the CMB, including the first signif- 
icant detection of the contribution from the SZ secondary 
anisotropy. 

In this paper we expand upon the results in S09 and 
present an SZ-detection-signiflcance-limited catalog of 
galaxy clusters identifled in the 2008 SPT survey. Red- 
shifts for 21 of these objects have been obtained from 
follow-up optical imaging, the details of which are dis- 
cussed in a companion paper (High et al. 2010). Using 
simulated observations we characterize the SPT cluster 
selection function — the detectability of galaxy clusters 
in the survey as a function of mass and redshift — for 
the 2008 fields. A simulation-based mass scaling relation 
allows us to compare the catalog to theoretical predic- 
tions and place constraints on the normalization of the 
matter power spectrum on small scales, ag, and the dark 
energy equation of state parameter w. 

This paper is organized as follows: §2 discusses the ob- 
servations, including data reduction, mapmaking, filter- 
ing, cluster-finding, optical follow-up and cluster redshift 
estimation; §3 presents the resulting cluster catalog; §4 
provides a description of our estimate of the selection 
function; §5 investigates the sample in the context of 
our current cosmological understanding and derives pa- 
rameter constraints; we discuss limitations and possible 
contaminants in §6, and we close with a discussion in §7. 

For our fiducial cosmology we assume a spatially flat 
ACDM model (parameterized by fifo/i^, Qch^, Hq, rig, r, 
and ^002) with parameters consistent with the WMAP 
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5-year ACDM best-fit results (Dunkley et al. 2009)^, 
namely Om = 0.264, = 0.044, h = 0.71, cts = 0.80. 
All references to cluster mass refer to M200, the mass 
enclosed within a spherical region of mean overdensity 
200 X Pmean, where Pmean is the mean matter density on 
large scales at the redshift of the cluster. 

2. OBSERVATIONS, DATA REDUCTION, CLUSTER 
EXTRACTION, AND OPTICAL FOLLOW-UP 

2.1. Observations 

The results presented in this work arc based on ob- 
servations performed by the SPT in 2008. Carlstrom 
et al. (2009) and S09 describe the details of these obser- 
vations; we briefly summarize them here. Two fields were 
mapped to the nominal survey depth in 2008: one cen- 
tered at right ascension (R.A.) 5'^30'", decUnation (dccl.) 
-55° (J2000), hereafter the 5^ field; and one centered at 
R.A. 23'^30'^, decl. -55°, hereafter the 23*^ field. Results 
in this paper are based on roughly 1500 hours of observ- 
ing time split between the two fields. The areas mapped 
with near uniform coverage were 91 deg^ in the 5^ field 
and 105 deg^ in the 23*^ field. 

This work considers only the 150 GHz data from the 
uniformly covered portions of the 2008 fields. The de- 
tector noise for the 95 GHz detectors was very high for 
the 2008 season, and the 220 GHz observations were con- 
taminated by the atmosphere at large scales where they 
would be \iscful for removing CMB fluctiiations. Includ- 
ing these bands did not significantly improve the effi- 
ciency of cluster detections and they were not used in the 
analysis presented here. The final depth of the 150 GHz 
maps of the two fields is very similar, with the white 
noise level in each map equal to 18 /xK-arcmin. 

The two fields were observed using slightly different 
scan strategies. For the 5^ field, the telescope was swept 
in azimuth at a constant velocity (~0.25°/s on the sky 
at the field center) across the entire field then stepped 
in elevation, with this pattern continuing until the whole 
field was covered. The 23^ field was observed using a 
similar strategy, except that the azimuth scans covered 
only one half of the field at any one time, switching halves 
each time one was completed. One consequence of this 
observing strategy was that a narrow strip in the middle 
of the 23^ field received twice as much detector time as 
the rest of the map. The effect of this strip on our catalog 
is minimal and is discussed in §4.6. 

A single observation of either field lasted ^ 2 hours. 
Between individual observations, several short calibra- 
tion measurements were performed, including measure- 
ments of a chopped thermal source, 2 degree elevation 
nods, and scans across the galactic HII regions RCW38 
and MAT5a. This series of regular calibration mea- 
surements was used to identify detectors with good per- 
formance, assess relative detector gains, monitor atmo- 
spheric opacity and beam parameters, and model point- 
ing variations. 

2.2. Data processing and mapmaking 

The data reduction pipeline applied to SPT data in 
this work is very similar to that used in S09. Broadly, 

^ These parameters arc sufficiently similar to the WMAP 7-ycar 
preferred cosmology (Larson et al. 2010) that a re-analysis based 
on that newer work is not warranted. 



the pipeline for both fields consists of filtering the time- 
ordered data from each individual detector, reconstruct- 
ing the pointing for each detector, and combining data 
from all detectors in a given observing band into a map 
by simple inverse-variance-weighted binning and averag- 
ing. 

The small differences between the data reduction used 
in this work and that of S09 arc: 

• In S09, a 19th-order polynomial was fit and re- 
moved from each detector's timestream on each 
scan across the field. Samples in the timestream 
which mapped to positions on the sky near bright 
point sources were excluded from the fit. A simi- 
lar subtraction was performed here, except that a 
first-order polynomial was removed, supplemented 
by sines and cosines (Fourier modes). Frequen- 
cies for the Fourier modes were evenly spaced from 
0.025 Hz to 0.25 Hz. This acts approximately as a 
high-pass filter in the R.A. direction with a char- 
acteristic scale of 1° on the sky. 

• In S09, a mean across functioning detectors was 

calculated at each snapshot in time and subtracted 
from each sample. Here, both a mean and a slope 
across the two-dimensional array were calculated at 
each time and subtracted. This acts as a roughly 
isotropic spatial high-pass filter, with a character- 
istic scale of ^ 0.5°. 

• As in Vieira et al. (2010), a small pointing cor- 
rection (^ 5" on the sky) was applied, based on 
comparisons of radio source positions derived from 
SPT maps and positions of those sources in the 
AT20G catalog (Murphy et al. 2010). 

The relative and absolute calibrations of detector re- 
sponse were performed as in LIO. The relative gains of 
the detectors and their gain variations over time were 
estimated using measurements of their response to a 
chopped thermal source. These relative calibrations were 
then tied to an absolute scale through direct comparison 
of WMAP 5-year maps (Hinshaw et al. 2009) to dedi- 
cated large-area SPT scans. This calibration is discussed 
in detail in LIO, and is accurate to 3.6% in the 150 GHz 
data. 

2.3. Cluster Extraction 

The cluster extraction procedure used in this work for 
both fields is identical to the procedure used in S09, 
where more details can be found. 

The SPT maps were filtered to optimize detection of 
objects with morphologies similar to the SZ signatures 
expected from galaxy clusters, through the application 
of spatial matched filters (Haehnelt & Tegmark 1996; 
Herranz et al. 2002a, b; Melin et al. 2006). In the spatial 
Fourier domain, the map was multiplied by 

^(fc.,fc,) = B{k.,ky)Si\k\) 

B{kx, ky)"^ Nastro{\k\) -\- Nnoise{kxj ky) 

where V is the matched filter, B is the response of the 

SPT instrument after timestream processing to signals 
on the sky, S is the assumed source template, and the 
noise power has been broken into astrophysical (Nastro) 
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and noise {Nnoise) components. For the source template, 
a projected spherical /3-model, with /3 fixed to 1, was 
used: 

AT = ATo(l + ^V^c)"^ 
whore the normalization ATq and the core radius Oc are 
free parameters. 

The noise power spectrum Nnoise includes contribu- 
tions from atmospheric and instrumental noise, while 
Nastro includes power from primary and lensed CMB 
fluctuations, an SZ backgroimd, and point sources. The 
atmospheric and instrumental noise were estimated from 
jackknife maps as in S09, the CMB power spectrum was 
updated to the lensed CMB spectrum from the WMAP5 
best-fit ACDM cosmology (Dunkley et al. 2009), the SZ 
background level was assumed to be flat in + 1)C^ 
with the amplitude taken from LIO, and the point source 
power was assumed to be flat in Ci at the level given in 
Hah et al. (2010). 

To avoid spurious decrements from the wings of bright 
point sources, all positive sources above a given flux 
(roughly 7 m Jy, or 5 cr in a version of the map filtered 
to optimize point-source signal-to-noise) were masked 
to a radius of 4' before the matched filter was applied. 
Roughly 150 sources were masked in each field, of which 
90-95% arc radio sources. The final sky areas considered 
after source masking were 82.4 and 95.1 deg^ for the 5^ 
and 23^ fields respectively. 

The maps were filtered for twelve different cluster 
scales, constructed using source templates with core radii 
9c evenly spaced between 0.25' to 3.0'. Each filtered map 
Mjj, where i refers to the filter scale and j to the field, 
was then divided into strips corresponding to distinct 
90' ranges in elevation. The noise was estimated inde- 
pendently in each strip in order to account for the weak 
elevation dependence of the survey depth. The noise in 
the fc*'' strip of map Mjj, a^/c, was estimated as the stan- 
dard deviation of the map within that strip. 

Signal-to-noise maps Mij were then constructed by 
dividing each strip k in map Mij by ajj^. SZ cluster 

decrements were identified in each map Mij by a simple 

(negative) peak dc;t(x;tion algorithm similar to SExtrac- 
tor (Bertin & Arnouts 1996). Tlic highest signal-to-noise 
value associated with a decrement, across all filter scales, 
was defined as ^, and taken as the significance of a detec- 
tion. Candidate clusters were identified in the data down 
to ^ of 3.5, though this work considers only the subset 
with ^> 5. 

These detection significances arc robust against choice 
of source template: the use of Nagai (Nagai et al. 2007), 
Arnaud (Arnaud et al. 2010) or Gaussian templates in 
place of ^—models was found to be free of bias and to 
introduce negligible (~ 2%) scatter on recovered ^. 

2.4. Optical Imaging and Spectroscopy 

Optical imaging was used for confirmation of candi- 
dates, for photometric redshift estimation, and for clus- 
ter richness characterization. A detailed description of 
the coordinated optical effort is presented in High et al. 
(2010) and is summarized here. 

The and 23^ fields were selected in part for over- 
lap with the Blanco Cosmology Survey (BCS, see Ngeow 
et al. 2009), which consists of deep optical images us- 
ing g, r, i and z filters. These images, obtained from 



the Blanco 4m telescope at CTIO with the MOSAIC- 
II imager in 2005 2007, were used when available. The 
co-added BCS images have 5 a galaxy detection thresh- 
olds of 24.75, 24.65, 24.35 and 23.05 magnitude in griz, 
respectively. 

For clusters that fell outside the BCS coverage region, 
as well as for 5 that fell within, and for the unconfirmed 
candidate, images were obtained using the twin 6.5 m 
Magellan telescopes at Las Campanas, Chile. The imag- 
ing data on Magellan were obtained by taking succes- 
sively deeper images until a detection of the early-type 
cluster galaxies was achieved, complete to between L* 
and 0.4L* in at least one band. The Magellan images 
were obtained under a variety of conditions, and the Stel- 
lar Locus Regression technique (High et al. 2009) was 
used to obtain precise colors and magnitudes. 

Spectroscopic data were obtained for a subset of the 
sample using the Low Dispersion Survey Spectrograph 
on Magellan (LDSS3, see Osip et al. 2008) in longslit 
mode. Typical exposures were 20-60 minutes, with slit 
orientations that contained the brightest cluster galaxy 
(BCG) and as many additional red cluster members as 
possible. 

Photometric redshifts were estimated using standard 
red sequence techniques and verified using the spectro- 
scopic subsample. A red sequence model was derived 
from the work of Bruzual & Chariot (2003), and local 
overdensities of red galaxies were searched for near the 
cluster candidate positions. By comparing the resulting 
photo-z's to spectroscopic redshifts within the subsample 
of 10 clusters for which spectroscopic data are available, 
HIO estimates photo- z uncertainty cr^ as given in Table 
1, of order 3% in (1 + z). 

Completeness in red sequence cluster finding was esti- 
mated from mock catalogs. At BCS depths, galaxy clus- 
ter completeness for the masses relevant for this sample 
is nearly unity up to a redshift of one, above which the 
completeness falls rapidly, reaching 50% at about redshift 
1.2, and 0% at redshift 1.25. At depths about a magni- 
tude brighter (corresponding to the depth of the Magel- 
lan observations of the unconfirmed candidate, §3.1), the 
completeness deviates from unity at redshift 0.8. 

3. CATALOG 

The resulting catalog of galaxy clusters, complete for 
^> 5, is presented in Table 1. Simulations (§4.1) sug- 
gest that this catalog should be highly complete above 
limiting mass and redshift thresholds, with relatively low 
contamination. 

A total of 22 candidates were identified, for which op- 
tical follow-up confirmed and obtained redshift informa- 
tion on all but one. Three clusters were previously known 
from X-ray and optical surveys, three were previously re- 
ported from this survey by S09, three were first identi- 
fied in a recent analysis of BCS data by Menanteau et al. 
(2010), and the remainder are new discoveries. Detailed 
comparisons of the SPT and Menanteau et al. (2010) 
cluster catalogs and selection will be the subject of fu- 
ture work. 

Thumbnail images of the signal-to-noise maps, M, at 
the preferred filter scale for each cluster are provided in 

Appendix A, Figure Al. Signal-to- noise as a function 
of filter scale for each cluster is shown in Appendix A, 
Figure A2. 
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TABLE 1 

The SPT Cluster Catalog for the 2008 Observing Season 
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Note. — Recall that ^ is the maximum signal-to-noise obtained over the set of filter scales for each cluster. The 6c (given 
in arcminutes) refer to the preferred filter scale, that at which ^ was found. Cluster positions in R.A. and decl. are given in 

degrees, and refer to the center of SZ brightness in map filtered at the preferred scale, calculated as the mean position of all 
pixels associated with the detection, weighted by their SZ brightness. The four rightmost columns refer to optical follow-up 
observations, giving the measured photometric redshift measurements and uncertainties of the optical counterpart, spectroscopic 
redshifts where available, and the source (BCS or Magellan) of the follow-up data. 

tClustcrs identified in S09. The cluster names have been updated (clusters were identified as SPT-CL 0509-5342, SPT-CL 0517-5430, 
SPT-CL 0528-5300, and SPT-CL 0547-5345 in S09) in response to an lAU naming convention request, an improved pointing model, and 
tlic updated data processing. 

tClustcrs within 2' of RASS sources (RASS-FSC, RASS-BSC; Voges et al. 2000, 1999). 

^SCSO J051145-515430 (Menanteau et al. 2010). 

''Abell S0520 (Abell et al. 1989), RXCJ0516.6-5430 (Bohringer et al. 2004), SCSO J051637-543001 (Menanteau et al. 2010). 
=SCSO J052113-510418 (Menanteau et al. 2010). 
<*SCSO J052803-525945 (Menanteau et al. 2010). 

''Abell S0552 (Abell et al. 1989) is in the foreground, 5' away at z=0.09 (this redshift not previously measured), 
f Abell 3950 (Abell et al. 1989). Spectroscopic redshift from Jones et al. (2005). 
^Abell S1079 (Abell et al. 1989, redshift shown not previously measured) . 
''SCSO J233227-535827 (Menanteau et al. 2010). 



Estimates of cluster masses are possible with the aid 
of a scaling relation (below, §4.8), and are discussed in 
Appendix D. 

3.1. Noteworthy Clusters 

SPT-CL 2259-5617— The SZ signal from this cluster is 
anomalously compact for such a low redshift object. The 
cosmological analysis (§-5) cixplicitly excludes all z < 0.3 
clusters, so this clustcir not used in parameter estimation. 

SPT-CL J2331-5051— This cluster appears to be one of 
a pair of clusters at comparable redshift, likely undergo- 
ing a merger. It will be discussed in detail in a future 
publication. The fainter partner is not included in this 
catalog as its significance = 4.81) falls below the de- 
tection threshold. 

SPT-CL J2332-5358— This cluster is coincident with a 
bright dusty point source which we identify in the 23h 
220 GHz data. Although the 150 GHz flux from this 



source could be removed with the aid of the 220 GHz 
map, a multi-frequency analysis is outside the scope of 
the present work. The impact of point sources on the 
resulting cluster catalog is discussed in §6.2. 

SPT-CL J2343-5521— No optical counterpart was found 
for this candidate. The field was imaged with both BCS 
and Magellan, and no cluster of galaxies was found to a 
5 a point source detection depth. The simulated optical 
completeness suggests that this candidate is either a false 
positive in the SPT catalog or a cluster at high redshift 
{z > 1.2). While the relatively high ^ = 5.74 indicates a 
^ 7% chance of a false detection in the SPT survey area 
(see discussion of contamination below, §4.7), the signal- 
to-noise of this detection exhibits peculiar behavior with 
9c (see Figure A2), preferring significantly larger scales 
than any other candidate, consistent with a CMB decre- 
ment. Further multi- wavelength follow-up observations 
are underway on this candidate, and preliminary results 
indicate it is likely a false detection. 
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3.2. Recovering integrated SZ flux 

The optimal filter described in §2.3 provides an esti- 
mate of the /3-model normalization, ATq, and core size 
for each cluster, based on the filter scale 9c at which the 
significance ^ was maximized. Assuming prior knowledge 
of the ratio 02oo/^c (where 6*200 is the angle subtend- 
ing the physical radius i?200 at the redshift of the clus- 
ter), one can integrate the /3-profile to obtain an estimate 
of the integrated SZ flux, Y. Basic physical arguments 
and hydrodynamical simulations of clusters have demon- 
strated F to be a tight (low intrinsic scatter) proxy for 
cluster mass (Barbosa et al. 1996; Holder & Carlstrom 
2001; Motl et al. 2005; Nagai et al. 2007; Stanek et al. 
2010). 

In single-frequency SZ surveys, the primary CMB tem- 
perature anisotropics provide a source of astrophysical 
contamination that greatly inhibits an accurate measure 
of (Melin et al. 2006). The modes at which the primary 
CMB dominates must be filtered out from the map, sig- 
nificantly reducing the range of angular scales that can be 
used by the optimal filter to constrain Oc- This range is 
already limited by the ^ 1' instrument beam, which only 
resolves 9c for the larger clusters. Any integrated quan- 
tity will thus be poorly measured. Melin et al. (2006) 
demonstrated that if the value of 9,, can be provided by 
external observations, e.g.. X-ray, Y can be accurately 
measured. 

The inability to constrain 9c can be seen in Figure A2, 
where the highest signal-to-noise associated with a peak 
is plotted against each filter scale. For several clusters 
(for example, SPT-CL J0516-5430, SPT-CL J0551-5709, 
and SPT-CL J2332-5358), the peak in signal-to-noise as- 
sociated with a cluster is very broad in 9c- Because of 
this confusion, we do not report Y in this work. Instead, 
as described below (§4.8), we use detection significance 
as a proxy for mass. 

Multi-frequency surveys are not in principle subject to 
this limitation as the different frequencies can be com- 
bined to eliminate sources of noise that are correlated be- 
tween bands, thus increasing the range of angular scales 
available for constraining cluster profiles. 

4. SZ SELECTION FUNCTION 

In this section, we characterize the SPT cluster sample 
identified in Table 1. Specifically, we describe the SPT 
cluster selection function in terms of the catalog com- 
pleteness as a function of mass and redshift, and the con- 
tamination rate. This selection function was determined 
by applying the cluster detection algorithm described in 
§2.3 to a large number of simulated SPT observations. 
These simulations included the dominant astrophysical 
components (primary and lensed CMB, cluster thermal 
SZ, and two families of point sources), accounted for the 
effects of the SPT instrument and data processing (the 
"transfer function" ) , and contained realistic atmospheric 
and detector noise. 

4.1. Simulated thermal SZ Cluster Maps 

Simulated SZ maps were generated using the method 
of Shaw et al. (2009), where a detailed description of the 
procedure can be formd. In brief, the semi-analytic gas 
model of Bode et al. (2007) was applied to halos identified 
in the output of a large dark matter lightcone simulation. 



The cosmological parameters for this simulation were 
chosen to be consistent with those measured from the 
WMAP 5-year data combined with large-scale structure 
observations (Dunkley et al. 2009), namely Qm = 0.264, 
fib = 0.044, and as = 0.8. The simulated volume was a 
periodic box of size 1 Gpc/h. The matter distribution in 
421 time slices was arranged into a lightcone covering a 
single octant of the sky from < z < 3. 

Dark matter halos were identified and gas distributions 
were calculated for each halo using the semi-analytic 
model of Bode et al. (2007). This model assumes that 
intra-clustcr gas resides in hydrostatic equilibrium in the 
gravitational potential of the host dark matter halo with 
a polytropic equation of state. As discussed in Bode 
et al. (2007), the most important free parameter is the en- 
ergy input into the cluster gas via non-thermal feedback 
processes, such as supernovae and outflows from active 
galactic nuclei ( AGN) . This is set through the parameter 
6/ such that the feedback energy is Ef — e/M^.c^, where 
is the total stellar mass in the cluster. Bode et al. 
(2007) calibrate ej by comparing the model against ob- 
served X-ray scaling relations for low redshift (z < 0.25) 
group and cluster mass objects. We note that the redshift 
range in which the model has been calibrated and that 
encompassed by the cluster sample presented here barely 
overlap; comparison of the model to the SPT sample (as, 
for example, in §5) thus provides a test of the predicted 
cluster and SZ signal evolution at high redshift. 

For our fiducial model we adopt e/ = 5 x 10~^, however 
for comparison we also generate maps using the 'stan- 
dard' and 'star-formation only' versions of this model 
described in Bode et al. (2009). There are two prin- 
cipal differences between these models. First, the stel- 
lar mass fraction M*/Mgas is constant with total clus- 
ter mass in the fiducial model, but mass-dependent in 
the 'standard' and 'star-formation' model. Second, the 
amount of energy feedback is significantly lower in the 
'standard' than model than in the fiducial, and zero in 
the 'star- format ion' model. 

From the output of each model, a 2-d image of SZ in- 
tensity for each cluster with mass M > 5 x 10^^ Mq/i"^ 
was produced by summing up the electron pressure along 
the line of sight. SZ cluster sky maps were constructed by 
projecting down the lightcone. summing up the contribu- 
tion of all the clusters along the line of sight. Individual 
SZ sky maps were 10 x 10 degrees in size, resulting in a 
total of 40 independent maps. For each map, the mass, 
redshift, and position of each chistcr was recorded. 

From SPT pointed observations of X-ray-selected clus- 
ters, Plagge et al. (2010) have demonstrated that cluster 
radial SZ profiles match the form of the "universal" elec- 
tron pressure profile measured by Arnaud et al. (2010) 
from X-ray observations of massive, low-rcdshift clusters. 
To complement the set of maps generated using the semi- 
analytic gas model, SZ sky maps were generated in which 
the projected form of the Arnaud et al. (2010) pressure 
profile was used to generate the individual cluster SZ 
signals. 

4.2. Point Source Model 

At 150 GHz and at the flux levels of interest to this 

analysis (^ 1 to ^ 10 mJy), the extragalactic source 
population is expected to be primarily composed of two 
broad classes: sources dominated by thermal emission 
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from dust heated by a burst of star formation, and 
sources dominated by synchrotron emission from AGN. 
We refer to these two famihes as "dusty sources" and 
"radio sources" and include models for both in our sim- 
ulated observations. 

For dusty sources, the source count model of Negrello 
et al. (2007) at 350 GHz was used. These counts arc 
based on the physical model of Granato ct al. (2004) for 
high-redshift SCUBA-like sources and on a more phe- 
nomenological approach for late-type galaxies (starburst 
plus normal spirals). Source counts were estimated at 
150 GHz by assuming a scaling for the flux densities of 
Si, oc v" , with a = 3 for high-redshift protospheroidal 
galaxies and a = 2 for late- type galaxies. For radio 
sources, the de Zotti et al. (2005) model for counts at 
150 GHz was used. This model is consistent with the 
measurements of Vicira et al. (2010) for the radio source 
population at 150 GHz. 

Realizations of source populations were generated by 
sampling from Poisson distributions for each population 
in bins with fluxes from 0.01 mJy to 1000 mJy. Sources 
were distributed in a random way across the map. Corre- 
lations between sources or with galaxy clusters were not 
modeled, and we discuss this potential contamination in 
§6.2. 

4.3. CMB Realizations 

Simulated CMB anisotropies were produced by gener- 
ating sky realizations based on the power in the gravita- 
tionally lensed WMAP 5-year ACDM CMB power spec- 
trum. Non-gaussianity in the lensed power was not mod- 
eled. 

4.4. Transfer Function 

The effects of the transfer function, i.e., the effects 
of the instrument beam and the data processing on 
sky signal, were emulated by producing synthetic SPT 
timestrcams from simulated skies sampled using the same 
scans employed in the observations. The sky signal 
was convolved with the measured SPT 150 GHz beam, 
timestream samples were convolved with detector time 
constants, and the SPT data processing (§2.2) was per- 
formed on the simulated timestreams to produce maps. 

Full emulation of the transfer function is a computa- 
tionally intensive process; to make a large number of sim- 
ulated observations, the transfer function was modeled 
as a 2D Fourier filter. The accuracy of this approxima- 
tion was measured by comparing recovered ^ of simulated 
clusters in skies passed through the full transfer function 
against the ^ of the same clusters when the transfer func- 
tion was approximated as a Fourier filter applied to the 
map. Systematic differences were found to be less than 
1% and on an object-by-object basis the two methods 
produced measured ^ that agreed to better than 3%. 

4.5. Instrum,ental & Atmospheric Noise 

Noise maps were created from SPT data by subtracting 
one half of each observation from the other half. Within 
each observation, one direction (azimuth either increas- 
ing or decreasing) was chosen at random, and all data 
when the telescope was moving in that direction were 
multiplied by -1. The data were then processed and com- 
bined as usual to produce a "jackknife" map which con- 



tained the full noise properties of the final field map, but 
with all sky signal removed. 

4.6. The 23^ Deep Strip 

Due to the observing strategy employed on the 23*^ 
field, a ~ 1.5° strip in the middle of that map contains 
significantly lower atmospheric and instrumental noise 
than the rest of the map. The jackknife noise maps 
(§4.5) used in simulated observations naturally include 
this deep strip, so any effects due to this feature are 
taken into account in the simulation-based estimation of 
the average selection function (§4.7) and scaling relation 
(§4.8) across the whole survey region. The cosmological 
analysis (§5) uses these averaged quantities; simulated 
observations performed with and without a deep strip 
demonstrated that any bias or additional scatter from 
using the averaged quantities is negligible compared to 
the statistical errors. 

4.7. Completeness and Contamination 

Forty realizations of the 2008 SPT survey (two fields 
each) were simulated, from which clusters were extracted 
and matched against input catalogs. 

Figure 1 shows the completeness of the simulated SPT 
sample, the fraction of clusters in simulated SPT maps 
that were detected with ^ > 5, as a function of mass 
and redshift. The exact shape and location of the curves 
in this figure depend on the detailed modeling of intra- 
cluster physics, which remain uncertain. The increase in 
SZ brightness (and cluster detectability) with increasing 
redshift at fixed mass is due to the increased density and 
temperature of high redshift clusters, and is in keeping 
with self-similar evolution. At low redshifts {z < 0.3), 
CMB confusion suppresses cluster detection significances 
and drives a strong low-redshift evolution in the selection 
function. These completeness curves were not used in 
the cosmological analysis (§5), where uncertainties on the 
mass scaling relation (§4.8) account for uncertainties in 
the modeling of intra-cluster physics. 

The SZ sky was removed from simulations to esti- 
mate the rate of false positives in the SPT sample. Fig- 
ure 2 shows this contamination rate as a function of 
lower ^ threshold, averaged across the survey area. A 
^ > 5 threshold leads to approximately one false detec- 
tion within the survey area. 

To test for biases introduced by an SZ background 
composed of low-mass systems, a simulation was run in- 
cluding only SZ sources well below the SPT threshold, 
with masses M < 10^* MQh~^. This background was 
found to have negligible effect on the detection rate as 
compared to the SZ-£ree false detection simulation. 

4.8. Mass Scaling Relation 

As discussed in §3.2, the integrated SZ flux Y is poorly 
estimated in this analysis and so is not used as a mass 
proxy. However, the noise aijk measured in each eleva- 
tion strip is relatively even across the SPT maps, so it 
is possible to work in the native space of the SPT selec- 
tion function and use detection significance ^ as proxy 
for mass. Additional uncertainty and bias introduced by 
use of such a relation (in place of, for example, a F-based 
scaling relation) are small compared to the Poisson noise 
of the sample and the uncertainties in modeling intra- 
cluster physics. 
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Fig. 1. — Simulated catalog completeness as a function of 
mass and redshift for a significance cut of g > 5. The contours 
show lines of constant completeness. From left to right, the lines 
represent 30, 50, 80 and 99% completeness. The temperature and 
density of clusters at a given mass tends to increase with redshift, 
leading to the increased SZ flux and improved detectability of 
high-redshift clusters. The strong evolution below z ~ 0.3 arises 
from reduced ^ on nearby clusters due to CMB confusion. Note 
that these contours are based on the fiducial simulations used 
in this work. Uncertainties in modeling (discussed in §4.8) can 
shift the position and shape of these contours coherently but 
significantly (of order 30% in mass). 




^ Lower Cutoff 

Fig. 2. — Simulated false detection rate, averaged across the 
survey area. The left axis shows the number density of false 
detections above a given ^; the right axis shows the equivalent 
number of false detections within the combined S*" and 23*^ survey 
fields. The dotted lines show the ^ > 5 threshold applied to the 
catalog, and the false detection rate at that threshold, ~ 1.2 
across the full survey area. 

The steepness of the cluster mass function in the pres- 
ence of noise will result in a number of detections that 
have boosted significance. Explicitly, ^ is a biased es- 
timator for (^), the average detection significance of a 
given cluster across many noise realizations. An ad- 
ditional bias on ^ comes from the choice to maximize 
signal-to- noise across three free parameters, R.A., decl. 
and dc- These biases make the relation between ^ and 
mass complex and difficult to characterize. 

In order to produce a mass scaling relation with a sim- 
ple form, the unbiased significance C is introduced. It is 
defined as the average detection signal-to-noise of a sim- 
ulated cluster, measured across many noise realizations, 



evaluated at the preferred position and filter scale of that 
cluster as determined by fitting the cluster in the absence 
of noise. 

Relating C, and ^ is a two step process. The expected 
relation between ^ and (^) is derived and compared to 
simulated observations in Appendix B, and found to be 
C = \l (^)^ ~ 3. Given a known (^), the expected distri- 
bution in ^ is derived by convolution with a Gaussian of 
unit width. The relation between C, and (^) is taken to be 
exact, and was verified through simulations to introduce 
negligible additional scatter; i.e., the scatter in the ^-^ 
relation is the same as the scatter in the (^)-^ relation, 
namely a Gaussian of unit width. 

The scaling between ^ and M is assumed to take the 
form of power-law relations with both mass and redshift: 

^^^(sx 1014Mq/i-i) (t6^) ' 

parameterized by the normalization A, the slope B, and 
the redshift evolution C. Appendix C presents a physical 
argument for the form of this relation, along with the 
expected ranges in which the values of the parameters B 
and C are expected to reside based on self-similar scaling 
arguments. 

Values for the parameters A, B, and C were deter- 
mined by fitting Eq. 1 to a catalog of C > 1 clus- 
ters detected in simulated maps, using clusters with 
mass M > 2 X 10^^ Mq/i"-'^ and in the redshift range 
0.3 < 2: < 1.2. This redshift range was chosen to match 
the SPT sample, while the mass limit was chosen to be 
as low as possible without the sample being significantly 
cut off by the C > 1 threshold. The best fit was de- 
fined as the combination of parameters that minimized 
the intrinsic fractional scatter around the mean relation. 

Figure 3 shows the best-fit scaling relation obtained 
for our fiducial simulated SZ maps, where A = 6.01, 
B = 1.31, and C = 1.6. The intrinsic scatter was mea- 
sured to be 21% (0.21 in ln(C)) and the relation was found 
to adhere to a power-law well below the limiting mass 
threshold. Over the three gas model realizations (§4.1), 
the best fit value of A, B, and the intrinsic scatter were 
all found to vary by less than 10%, while the values of C 
predicted by the 'standard' and 'star-formation' models 
drop to ~ 1.2. For maps generated using the electron 
pressure profile of Arnaud et al. (2010), best-fit values of 
A = 6.89, B = 1.38, C = 0.6 were found, with a 19% in- 
trinsic scatter. The values of A and B remain within 15% 
of the fiducial model, although C is significantly lower. 
Arnaud et al. (2010) measured the pressure profile using 
a low-redshift {z < 0.2) cluster sample and assume that 
the profile normalization will evolve in a self-similar fash- 
ion. The mass dependence of their pressure profile was 
determined using cluster mass estimates derived from the 
equation of hydrostatic equilibrium; simulations suggest 
that this method may underestimate the true mass by 
10 - 20% (Rasia et al. 2004; Meneghetti et al. 2010; Lau 
et al. 2009). We do not take this effect into account in 
our simulations - doing so would reduce the value of A by 
approximately 10%. The Bode et al. (2009) gas model is 
calibrated against X-ray scaling relations measured from 
low-redshift cluster samples (Vikhlinin et al. 2006; Sun 
et al. 2009), but assumes an evolving stellar- mass frac- 
tion which may drive the stronger redshift evolution. 
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Fig. 3. — Mass-significance relation plotted over clusters 
identified in simulated maps. The relation was fit to points with 
iU > 2 X 10^^ Mrj-Ji^^ , shown by the dotted line, and across a 
redshift range 0.3 < z < 1.2. Simulated clusters outside this 
redsliift range are not included in this plot. The approximate 
lower mass threshold of the high-redshift end of the SPT sample 
(M = 4 X lO^"' MQh~^) is shown by the dashed line. 

Based on these simulations, priors on the scahng rela- 
tion parameters {A, B, C, scatter) were adopted, with 
conservative la Gaussian uncertainties of (30%, 20%, 
50%, 20%) about mean values measured from the fiducial 
simulation model. These large uncertainties in scaling re- 
lation parameters arc the dominant source of uncertainty 
in the cosmological analysis (§5) and mass estimation 
(Appendix C). Furthermore, although the weakest prior 
is on the redshift evolution, C, it is the uncertainty on 
the amplitude A that dominates the error budget on the 
measurement of cig (sec §5.4). 

It should be noted that at low redshift (z < 0.3), such 
a power-law scaling relation fails to fully capture the be- 
havior of the CMB-confused selection function. The cos- 
mological analysis below therefore excludes this region 
during likelihood calculation. The mass estimates pre- 
sented in Appendix D may be biased low for low-redshift 
objects, although this effect is expected to be small com- 
pared to existing systematic errors. 

5. COSMOLOGICAL ANALYSIS 

The 2008 SPT cluster catalog is an SZ-detection- 
significance-limited catalog. Simulated maps were used 
to calibrate the statistics of the relation between cluster 
mass and detection significance, as well as the impact of 
noise-bias and selection effects. This relation was com- 
bined with theoretical mass functions to construct esti- 
mates of the number density of galaxy clusters as a func- 
tion of the significance $, and redshift, to be compared 
to the SPT catalog. Cosmological information from the 
SPT cluster catalog was combined with information from 
existing data sets, providing improved parameter con- 
straints. 

5.1. Cosmological Likelihood Evaluation 

Evaluation of cosmological models in the context of the 
SPT catalog requires a theoretical model that is capable 
of predicting the number density of dark matter halos 
as a function of both redshift and input cosmology. For 
a given set of cosmological parameters, the simulation- 
based mass function of Tinker et al. (2008) was used 



in conjunction with matter power spectra computed by 
CAME (Lewis et al. 2000) to construct a grid of clus- 
ter number densities in the native ^-z space of the SPT 
catalog: 

• A 2D grid of the number of clusters as a function of 

redshift and mass was constructed by multiplying 
the Tinker et al. (2008) mass function by the co- 
moving volume element. The gridding was set to be 
very fine in both mass and redshift, with Az = 0.01 
and the mass binning set so that = 0.0025 (see 
below). The grids were constructed to extend be- 
yond the sensitivity range of SPT, 0.1 < z < 2.6 
and 1.8 < C < 23. Extending the upper limits was 
found not to impact cosmological results, as pre- 
dicted number counts have dropped to negligible 
levels above those thresholds. 

• The parameterized scaling relation (§4.8) was used 
to convert the mass for each bin to unbiased signif- 
icance C for assumed values oi A,B and C. 

• This grid of number counts (in Q—z space) was con- 
volved with a Gaussian in ln((^) with width set by 
the assumed intrinsic scatter in the scaling relation 
(21% in the fiducial relation). 

• The unbiased significance Q of each bin was con- 
verted to an ensemble-averaged significance {£). 

• This grid was convolved with a unit-width Gaus- 
sian in ^ to account for noise, with the resulting 
grid in the native SPT catalog space, ^—z. 

• Each row (fixed ^) of the ^— z grid which contained 
a cluster was convolved with a Gaussian with width 
set as the redshift uncertainty for that cluster. 
Photometric redshift uncertainties are given in Ta- 
ble 1, and are described briefly in §2.4 and in detail 
in High et al. (2009); spectroscopic redshifts were 
taken to be exact. 

• A hard cut in ^ was applied, corresponding to the 
catalog selection threshold of ^ > 5. 

• An additional cut was applied, requiring z > 0.3, 
to avoid low-redshift regions where the power-law 
scaling relation fails to capture the behavior of 

the CMB-confused selection function. This cut ex- 
cludes 3 low-redshift clusters from the cosmological 
analysis, leaving 18 clusters plus the unconfirmed 
candidate, whose treatment is described below. 

The likelihood ratio for the SPT catalog was then con- 
structed, as outlined in Cash (1979), using the Poisson 
probability. 



N 



where the product is across bins in ^— z space, N is the 
number of bins. Pi is the Poisson probability in bin i, 
and Ci and rn are the fractional expected and integer 
observed number counts for that bin, respectively. 

The unconfirmed candidate was accounted for by si- 
multaneously allowing it to either be at high redshift or 
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a false detection. Its contribution to the total likelihood 
was calculated as the union of the likelihoods for n = 
and n = 1 within a large z > 1.0 bin, the redshift range 
corresponding to where the optical completeness for this 
candidate's follow-up deviates from unity. 

Ultimately, two sources of mass-observable scatter - 
the intrinsic scatter in the scaling relation, and the 1 a 
measurement noise - were included in this analysis, along 
with redshift errors and systematic uncertainties on scal- 
ing relation parameters. Other sources of bias and noise 
(such as point source contamination, §6.2, and the mass 
function normalization described below) are thought to 
be subdominant to these and were disregarded. 

While Tinker et al. (2008) claim a very small (< 5%) 
uncertainty in the mass function normalization, Stanek 
et al. (2010) have demonstrated that the inclusion of non- 
gravitational baryon physics in cosmological simulations 
can modify cluster masses in the range of the SPT sam- 
ple by ^ ±10% relative to gravity-only hydrodynamical 
simulations. The large 30% uncertainty on the ampli- 
tude A of the scaling relation effectively subsumes such 
uncertainties in the mass function normalization. 

This analysis does not account for the effects of sam- 
ple variance (Hu & Kravtsov 2003); for the mass and 
redshift range of the SPT sample, this is not expected 
to be a problem (Hu & Cohn 2006). The SPT survey 
fields span of order 100 h^^ Mpc, where the galaxy clus- 
ter correlation function would be expected to be a few 
percent or less (Bahcall et al. 2003; Estrada et al. 2009). 
This leads to clustering corrections to the uncertainty on 
number counts on the order of a few percent or less of 
the Poisson variance. 

5.2. Application to MCMC Chains 

The present SPT sample only meaningfully constrains 
a subset of cosmological parameters, so to explore the 
cosmological implications of the SPT catalog it is nec- 
essary to include information from other experiments. 
Existing analyses of other cosmological data in the form 
of Markov Chain Monte Carlo (MCMC) chains provide 
fully informative priors. These were importance sampled 
by weighting each set of cosmological parameters in the 
MCMC chain by the hkelihood of the SPT cluster catalog 
given that set of parameters. 

In this analysis, four MCMC chains were used to ex- 
plore constraints on parameters: the first two use only 
the 7-year data set from the WMAP experiment to ex- 
plore the standard spatially flat ACDM and wCDM cos- 
mologies, while the third explores wCDM while adding 
data from baryon acoustic oscillations (BAO) (Percival 
et al. 2010) and supernovae (SNe) (Hicken et al. 2009). 
These three chains were taken from the official WMAP 
analysis'^ (Komatsu et al. 2010). The fourth chain was 
computed by LIO and allows for a direct comparison with 
that work. It explores a spatially flat ACDM parameter 
space based on the "CMBall" data set: WMAP 5-year -f 
QUaD (Brown et al. 2009) -I- ACBAR (Reichardt et al. 
2009) -I- SPT (power spectrum measurements with Asz 
as a free parameter; LIO). 

Figure 4 shows the number density of clusters in the 
SPT catalog, plotted over theoretical predictions calcu- 
lated using the method described in §5.1, for 100 ran- 

^ Chains available at http://lambda.gsfc.nasa.gov 
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Fig. 4. — The SPT catalog, binned into 3 redshift bins 
(z=0.1-0.5, 0.5-0.9, 0.9-1.3), with number counts derived from 100 
randomly selected points in the WMAP7 wCDM MCMC chain 
overplotted. The SPT data are well covered by the chain and 
provide improved constraining power. The unconfirmed candidate 
is not included in this plot, and the binning is much coarser for 
display purposes than that used in the likelihood calculation (§5.1). 

dom positions in the WMAP7-only MCMC chain. The 
SPT data are adequately described by many cosmolog- 
ical models that are allowed by this data set, and the 
MCMC chains are well-sampled within the region of high 
probability. 

Uncertainties in the scaling relation parameters were 
accounted for by marginalizing over them: at each step 
in the chain, the likelihood was maximized across A, B, C 
and scatter, subject to the priors applied to each param- 
eter, using a Newton-Raphson method. The parameter 
values selected in this way at the highest likelihood point 
in each MCMC chain are given in Table CI. The fidu- 
cial values of B and the scatter appear consistent with 
those preferred by the chains, while the preferred values 
of the normalization A and redshift evolution C are both 
approximately 10% lower than their fiducial values. 

This weaker-than-fiducial redshift evolution could 
come from a variety of sources, and is consistent with 
other simulated models, e.g., the 'standard' and 'star- 
formation' models, see §4.8. Uncertainties in the redshift 
evolution are not a significant source of error in this anal- 
ysis: recovered parameter values and uncertainties (§5.3) 
are found to be insignificantly affected by widely varying 
priors on C. 

5.3. Cosmological Parameter Constraints 

The resulting constraints on ag, and w are given for 
all chains in Table 2. The parameter best constrained 
by the SPT cluster catalog is as- CMB power spectrum 
measurements alone have a large degeneracy between the 
dark energy equation of state, w, and fjg- Figure 5 shows 
this degeneracy, along with the added constraints from 
the SPT cluster catalog. Including the cluster results 
tightens the erg contours and leads to an improved con- 
straint on w. This is a growth-based determination of 
the dark energy equation of state, and is therefore com- 
plementary to dark energy measurements based on dis- 
tances, such as those based on SNe and BAO. 

When combined with the wCDM WMAP7 chain, the 
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Fig. 5. — Likelihood contour plot of w versus ag showing 1 a and 2 cr contours for several data sets. The left panel shows the constraints 
from WMAP7 alone (blue) and with the SPT cluster catalog included (red). The right panel shows show the full cosmological data set of 
WMAP7+SN+BAO (blue), and this plus the SPT catalog (red). The ability to constrain cosmological parameters is severely impacted 
by the uncertainties in the mass scaling relation, though some increase in precision is still evident. 



TABLE 2 

Cosmological Parameter, Constraints 



Chain 




^8 


w 


ACDM 


WMAP7 


0.801 ±0.030 


-1 


ACDM 


WMAP7+SPT 


0.791 ±0.027 


-1 


ACDM 


CMBall 


0.794 ±0.029 


-1 


ACDM 


CMBall+SPT 


0.788 ±0.026 


-1 


wCDM 


WMAP7 


0.832 ± 0.134 


-1.118 ±0.394 


wCDM 


WMAP7+SPT 


0.810 ±0.090 


-1.066 ±0.288 


wCDM 


WMAP7+BA0+SNe 


0.802 ±0.038 


-0.980 ±0.053 


wCDM 


WMAP7+BA0+SNe+SPT 


0.790 ±0.034 


-0.968 ± 0.049 



Note. — Mean values and symmetrized la range for erg and w, as found from each of the four data sets considered, shown 
with and without the weighting by hkehhoods derived from the SPT cluster catalog. The parameter best constrained by the 
SPT cluster catalog is ag. CMB power spectrum measurements alone have a large degeneracy between the dark energy equation 
of state, w, and o-g. Adding the SPT cluster catalog breaks this degeneracy and leads to an improved constraint on w. The 
SPT catalog has negligible effect on other parameters in these chains {fibh^, flch^, Ho, r and ris). 
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SPT data provide roughly a factor of 1.5 improvement 
in the precision of erg and w, finding 0.81 ± 0.09 and 
— 1.07 ± 0.29, respectively. Including data from BAO 
and SNe, these constraints tighten to as — 0.79 ± 0.03 
and w = -0.97 ±0.05. 

The dominant sources of uncertainty limiting these 
constraints are the Poisson error due to the relatively 
modest size of the current catalog and the uncertainty in 
the normalization A of the mass scaling relation. With 
weak-lensing- and X-ray-derived mass estimates of SPT 
clusters, along with an order of magnitude larger sample 
expected from the full survey, cosmological constraints 
from the SPT galaxy cluster survey will markedly im- 
prove. 

5.4. Amplitude of the SZ Effect 

The value of the normalization parameter A (which 
can be thought of as an "SZ amplitude" ) preferred by the 
likelihood analysis was found to be lower than the fidu- 
cial value, as shown in Figure 6. The prior assumed on 
this parameter is sufficiently large that it is not a highly 
significant shift; however, in light of the recent report by 
LIO of lower-than-expected SZ flux, it is worth address- 
ing. The SPT cluster catalog results are complementary 
to the the results of the power spectrum analysis, in that 
the majority of the SZ power at the angular scales probed 
by LIO comes from clusters below the mass threshold of 
the cluster catalog. 

Figure 6 shows that the amplitude A is strongly de- 
generate with (jg. The constraints provided by the SPT 
cluster catalog indicate either a value of as that is at 
the low end of the CMB-allowed distribution (or equiva- 
lently an erroneously high mass function normalization) , 
or an over-prediction of SZ flux by the fiducial simula- 
tions. If the fiducial amplitude is assumed, the best-fit as 
drops from the WMAP5+CMBaU value of 0.794 ± 0.029 
to 0.775 ± 0.015. This value is anomalously low com- 
pared to recent results (e.g. Vikhlinin et al. 2009; Mantz 
et al. 2010), and in slight tension with the results of the 
power spectrum analysis of LIO, where a still lower value 
of as = 0.746 ± 0.017 was obtained for similar simula- 
tion models.^ The SZ amplitude parameter used in LIO, 
Agz, is roughly analogous to A^ in the current notation. 
When including the expected contribution from homoge- 
neous reionization, LIO found Agz = 0.42 ± 0.21, in mild 
tension (at the ~ la level) with the marginalized value 
of (A/Afid)"^ = 0.79 ± 0.30 found in this analysis. 

The fiducial simulations in this work use the semi- 
analytic gas model of Bode et al. (2007, 2009), which 
is calibrated against low-redshift (z < 0.25) X-ray obser- 
vations but has not previously been compared to higher 
redshift systems. One interpretation of these results is 
that this model may over-predict the thermal electron 
pressure in high-redshift {z > 0.3) systems; this is not 
in conflict with the low-redshift calibration of the model 
and suggests a weaker redshift evolution in the SZ signal 
than predicted by the model. Alternately, a combination 
of mass function normalization and point source contam- 
ination could potentially account for the difference. 

The fiducial thermal SZ simulation model used in this paper 
predicts a power spectrum that is in very close agreement with the 
fiducial model of LIO, which was measured from the simulations of 
Sehgal et al. (2010). 
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Fig. 6. — Degeneracy between erg and SZ scaling relation 
amplitude A, plotted without prior (green) and with a 30% 
Gaussian prior (red) on A, for the ACDM WMAP5+CMBall 
MCMC chain. The Gaussian prior is shown (ztlcr) by the gray 
band, with the fiducial relation amplitude shown by the blue line. 
This figure is analogous to Fig. 9 of LIO, although that work dealt 
with SZ power, which is roughly proportional to the square of the 
amplitude being considered here. The prior is slightly higher than 
the preferred value; these results suggest that simulations may 
over-estimate the SZ flux in the high-mass, high-redshift systems 
contained in this catalog. 

6. SOURCES OF SYSTEMATIC UNCERTAINTIES 

There are several systematic effects that might affect 
the utility of the SPT cluster sample. For example, there 
remains large uncertainty in the mapping between detec- 
tion significance and cluster mass. It is also possible that 
strong correlations (or anti-correlations) between galaxy 
clusters and mm-bright point sources are significant. We 
address these issues in this section. 

6.1. Relation between SZ signal and Mass 

Theoretical arguments (Barbosa et al. 1996; Holder & 
Carlstrom 2001; Motl et al. 2005) suggest that the SZ flux 
of galaxy clusters is relatively well understood. However, 
there is very little high-precision empirical evidence to 
confirm these arguments, and there are physical mech- 
anisms that could lead to suppressed SZ fiux, such as 
non-thermal pressure support from turbulence (Lau et al. 
2009) or non-equilibrium between protons and electrons 
(Fox & Loeb 1997; Rudd & Nagai 2009). 

Cluster SZ mass proxies (such as Y and yo, the inte- 
grated SZ fiux and amplitude of the SZ decrement, re- 
spectively) depend linearly on the gas fraction and the 
gas temperature. There remain theoretical and obser- 
vational uncertainties in both of these quantities. Esti- 
mates of gas fractions for individual clusters can disagree 
by nearly 20% (e.g., Allen et al. 2008; Vikhlinin et al. 
2006), while theoretical and observed estimates of the 
mass-temperature relation currently agree at the level of 
10-20% (Nagai et al. 2007). Adding these in quadrature 
leads to uncertainties slightly below our assumed prior 
uncertainty of 30%. 

With the number counts as a function of mass, 
dN/dlnM, scaling as M'^ or for typical SPT clus- 
ters (Shaw et al. 2010a), a 10% offset in mass would lead 
to a 20-30% shift in the number of galaxy clusters. With 
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a catalog of 22 clusters, counting statistics lead to an un- 
certainty of at least 20%. Therefore, systematic offsets in 
the mass scale of order 10% will have a significant effect 
on cosmological constraints, and the current 30% prior 
on A will dominate Poisson errors. 

A follow-up campaign using optical and X-ray observa- 
tions will buttress our current theory/simulation-driven 
understanding of the SPT SZ-selected galaxy cluster cat- 
alog. 

6.2. Clusters Obscured by Point Sources 

The sky density of bright point sources at 150 GHz is 
low enough — on the order of 1 deg~^ (Vieira et al. 
2010) — that the probability of a galaxy cluster being 
missed due to a chance superposition with a bright source 
is negligible. However, sources associated with clusters 
will preferentially fill in cluster SZ decrements. Charac- 
terizing the contamination of cluster SZ measurements 
by member galaxies will be necessary to realize the full 
potential of the upcoming much larger SPT cluster cat- 
alog, but the systematic uncertainty predicted here and 
in the literature is well below the statistical precision of 
the current sample; it is disregarded in the current cos- 
mological analysis (§5). 

6.2.1. Dusty Source Contamination 

Star formation is expected to be suppressed in clus- 
ter environments (e.g., Hashimoto et al. 1998). Bai et al. 
(2007) measure the abundance of infrared- luminous star- 
forming galaxies in a massive (> IO^^Mq) cluster at 
z = 0.8 to be far lower relative to the field abundance 
than a simple mass scaling would predict: the cluster vol- 
ume that is hundreds of times ovcrdense in mass is only 
20 times ovcrdense in infrared luminosity. A sphere at 
z = 1 with a 1 Mpc radius and infrared luminosity that is 
20 times larger than the field would produce < 0.1 mJy of 
emission at 150 GHz, according to the sub- mm luminos- 
ity measurements of BLAST (Pascale et al. 2009) (Pascale 
et al, 2009). Even if the IR overdensity evolves strongly 
with mass and redshift, we can expect <C 1 mJy of con- 
tamination for the highest-redshift {z ~ 1) clusters at the 
SPT mass threshold. This corresponds to ^ 10% of the 
cluster SZ signal, which is far less than the uncertainty 
in the normalization of cluster masses presented in this 
work. 

Additionally, Keisler (2010, in prep.) measures the av- 
erage 100 /xm flux of cluster members from a sample of 
clusters at (z) — 0.2 and with masses similar to those se- 
lected by SPT and, after extrapolating to 150 GHz and 
allowing for strong redshift evolution in the infrared lumi- 
nosity function, constrains this contamination to be less 
than 10% of the cluster SZ signal. Again, this level of 
contamination is subdominant to the uncertainty in the 
normalization of cluster masses presented in this work. 

6.2.2. Gravitational Lensing 

Galaxy clusters can gravitationally lens sources located 
behind them. Because gravitational lensing conserves 
surface brightness, this process cannot alter the mean 
flux due to the background sources when averaged over 
many clusters. The background of sources is composed 
of both overdcnsities and underdensities, leading to both 
positive and negative fluctuations, relative to the mean, 
which will be gravitationally lensed. 



We do not explicitly account for this effect in this 
work. The unlenscd fluctuating background of sources at 
150 GHz is expected to be small (Hall et al. 2010) com- 
pared to both the experimental noise and intrinsic scatter 
on the mass scaling relation, and lensing only marginally 
increases the noise associated with these background 
sources (Lima et al. 2010). Within the context of the 
cosmological analysis, this additional noise term is ex- 
pected to be small compared to the intrinsic scatter on 
the mass scaling relation. 

6.2.3. Radio Source Contamination 

Galaxy clusters are known to host radio sources, but 
these correlated sources are not expected to be a major 
contaminant at 150 GHz. Calculations (Lin et al. 2009) 
and explicit simulations (Sehgal et al. 2010) demonstrate 
that, even taking into account the expected correlation 
between clusters and radio sources, these sources are not 
expected to significantly affect the SZ flux in more than 
1% percent of galaxy clusters above 2 x 10^** at a 
redshift of ^; ~ 0.5 (where "significantly" here means at 
the > 20% level). 

Simulations were also performed rising knowledge of 
the radio source popiilation at 150 GHz from Vieira ct al. 
(2010) and the cluster profiles that maximize the signif- 
icance for the SPT clusters presented here. Each profile 
between rcorc = 0.25' and rcoro = 1-5' (a range which en- 
compasses all of the optically confirmed clusters in Table 
1) was scaled so that the filtered version of that pro- 
file would result in a ^ = 5 detection in the 2008 SPT 
maps. Point sources of a given flux were then added 
at a given radius from the profile center. These point- 
source-contaminated profiles were then convolved with 
the transfer function, the matched filter was applied, 
and the resulting central value was compared to the cen- 
tral value of the filter-convolved, uncontaminated profile. 
Clusters were found to suffer a systematic = 1 reduc- 
tion in significance from a 2 mJy(5 mJy) source at 0.5'(1') 
from the profile center. This effect is nearly independent 
of core radius in the range of core radii probed. 

The Vieira et al. (2010) radio source counts at 150 GHz 
indicate roughly 1.5 per deg^ above 5mJy, while the de 
Zotti et al. (2005) 150 GHz model predicts roughly 3 
radio sources per deg^ above 2 mJy.^ If there were no 
correlation between clusters and radio sources, the clus- 
ters contained in the SPT catalog should have a 0.14% 
(0.03%) chance of incurring an error of > 1 from a 
> 5 mJy (2 — 5 mJy) source. 

Furthermore, using 30 GHz observations of a sample 
of clusters ranging from 0.14 < z < 1.0, Coble et al. 
(2007) find the probability of finding a radio source near a 
cluster to be ^-^^.2% (3.3l^;8) times the background rate 
when using a 0.5' (5') radius. From these results, it can 
be estimated that roughly 1% of SPT-detected clusters 
would suffer an error of A^ > 1 from radio source con- 
tamination. This is in very close agreement with the pre- 
dictions firom Lin et al. (2009) and Sehgal et al. (2010). 

7. DISCUSSION 

^ The Vieira et al. (2010) counts do not cover a low enough flux 
range to predict counts at 2 m,Jy, but the de Zotti ct al. (2005) 
model is consistent with the Vieira ct al. (2010) counts at all fluxes 
above 5 mJy, so counts from this model can confidently be extrap- 
olated down a factor of 2.5 in flux. 
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We have presented the first cosmologicahy significant 
SZ-selected galaxy cluster catalog, characterized the se- 
lection function, and performed a preliminary cosmologi- 
cal analysis to both demonstrate the general consistency 
of the catalog with current understanding of cosmology 
and provide improved constraints on cosmological pa- 
rameters. This is an important step toward exploiting 
the potential of SZ-sclcctcd galaxy clusters as a powerful 
cosmological tool. 

Using single-frequency data taken in 2008 with SPT, 
a total of 22 candidates were identified, all but one of 
which were optically confirmed as galaxy clusters (High 
et al. 2010). Of these 21 clusters, three were previously 
known from optical and/or X-ray surveys, three were new 
SPT detections reported in S09, three were first identified 
from BCS data by Menanteau et al. (2010), and 12 are 
new discoveries. 

Simulations were used to calibrate the selection func- 
tion of the survey and measure a scaling between SPT 
detection significance and mass. These simulations indi- 
cate that SZ detection significance traces mass with little 
(^ 20%) intrinsic scatter, making SZ surveys well suited 
to selecting mass-limited catalogs of galaxy clusters. 

As a demonstration of the constraining power of the 
survey, the SPT cluster catalog was used to refine es- 
timates of cosmological parameters, including the dark 
energy equation of state, w. and the normalization of 
the matter power spectrum on small scales, erg. Using 
wCDM MCMC chains derived from the WMAP 7-year 
data combined with the SPT cluster catalog, the best-fit 
values were u; = -1.07 ± 0.29 and as = 0.81 ± 0.09, a 
factor of roughly 1.5 improvement in precision compared 
to the WMAP7 constraints alone. When combined with 
other cosmological data sets (baryon acoustic oscillations 
and supernovae), the SPT cluster catalog improves pre- 
cision on these parameters by ~ 10%. 

These results can be compared to those of Vikhlinin 
et al. (2009) and Mantz et al. (2010), who performed 
a similar analysis using large samples of clusters drawn 
from X-ray surveys. The SPT results are less pre- 
cise: in combination with various cosmological data sets, 
Vikhlinin et al. (2009) find nearly 4 times tighter con- 
straints on cTg, while Mantz et al. (2010) arc more pre- 
cise by nearly a factor of 2. This is not surprising: both 
X-ray analyses had significantly larger cluster samples 
and smaller stated uncertainty in the mass scaling rela- 
tion. The weaker parameter constraints foimd from the 
SPT cluster catalog arc a direct result of uncertainties in 
the mass scaling relation, which derive from uncertainties 
modeling intra-cluster physics. 

The fiducial thermal SZ simulation model assumed 
here was shown to produce some tension between the 
analysis presented here and contemporary cosmological 
results. This may be explained by a variety of factors. 
The value of as may be lower than currently favored, 
or equivalently the normalization of the Tinker et al. 
(2008) mass function may be erroneously high. Alterna- 
tively, current simulations, while reproducing observed 
low-redshift X-ray observations, may over-estimate SZ 
flux in higher redshift systems, implying missing physics 
in the semi-analytic gas modeling (for example, a non- 
negligible amount of non-thermal pressure support at 
higher redshifts, Shaw et al. 2010b). The observed 
SZ signal could potentially be contaminated by an in- 



creasing incidence of point-source emission at high red- 
shift although the arguments presented in §6.2 suggest 
point-source contamination is unlikely to be wholly re- 
sponsible. LIO also found lower-than-anticipated power 
from the SPT measurement of the SZ power spectrum, 
consistent with many of these scenarios. A concerted, 
multi-wavelength program aimed at studying high red- 
shift clusters should help to resolve these issues. 

The SPT catalog presented here is based on less than 
1/3 of the current data and roughly 1/10 of the full SPT 
survey. The large multifrcquency SPT survey, combined 
with X-ray and/or weak lensing mass estimates of a sub- 
sample of SZ-selected galaxy clusters, should allow an 
order of magnitude improvement in the precision of as 
and w measurements. 
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APPENDIX 

A. THUMBNAIL IMAGES OF CLUSTER SZ DECREMENTS 




Fig. A1. — Thumbnail cutouts of filtered maps, 25' x 25', of the highest significance filter scale for each cluster, centered on the cluster 
position. The color scale and spatial extent of all maps here are the same, and the irregular shape of some images shows the edge of the 
masked (even coverage) map. The large blur in the upper right of SPT-CL J2331-5051 is a point source mask. 
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Fig. A2. — The measured significance across filter scales for all clusters. The maximum value of each curve gives ^ for that candidate. 
The significance axes have been allowed to float for each panel. Adjacent points are highly correlated, and as discussed in §3.2, the best-fit 
0c is poorly constrained. 
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B. UNBIASED SIGNIFICANCE C 

The freedom to maximize the SPT significance ^ across three parameters (R.A., decl., and 6c) in the presence of 

a noise field will tend to raise the amplitude of the observed peak. That is, the ensemble average of ^ across many 
noise realizations, (^), will be boosted by some amount as compared to the unbiased significance C, which is measured 
without these degrees freedom. 

We can consider the SPT significance ^ as analogous to a of the best-fit model relative to zero signal, with the 
signal to noise being comparable to \/x^. By allowing three degrees of freedom in a usual fit, we expect that the 
minimum will typically be smaller than the corresponding one for a fit with no degrees of freedom, i.e., at the true 
location using the true filter size. The typical difference is simply the number of fit parameters, in our case three: 

We test the relation on simulated data, adding 25 different realizations of each noise term (CMB, point sources, 
instrumental and atmospheric) to the same SZ sky, and explicitly measuring (^) for many simulated clusters in these 
maps. We next simulate and measure (. for the same patch of sky and compare the two measures of significance. The 
results of these simulations are shown in Figure Bl, and are consistent with the hypothesized relation with = 27 
for 23 degrees of freedom. 

We adopt this relation for the cosmological analysis and mass estimation, 

and find that any residual behavior not modeled by this relation is small enough to be well sub-dominant to uncertainties 
on the scaling relation amplitude A and Poisson noise within the catalog. Within the SPT catalog this is a small effect: 
for (^) = 5(10) objects, the difference between ( and (^) is 6%(1%). 
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Fig. Bl. — The quadratic difference between measured f and ($) (measured across 25 realizations of all noise terms), binned and plotted 
against a wide range of C- These data are consistent with {^)^ — = 3, with = 27 for 23 degrees of freedom. 
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C. MASS-SIGNIFICANCE SCALING PARAMETERS 
We have assumed that the scahng between C and M takes the form 



where chister mass M is defined in terms of a spherieal overdensity relative to the mean density at a given redshift 
z and the pivot points in this relation were chosen to roughly match the mean mass and redshift of clusters in the 
SPT sample. Estimates for B and C were obtained from extensive simulations; in this appendix we explore various 
limiting cases to better understand what values are expected. For the cosmological analysis we use priors informed by 
our simulations. 

The amplitude of the SZ effect through the center of a cluster (the central decrement) is expected to scale as (Holder 
& Carlstrom 1999) 

yoocM{l + zf , (C2) 

so if C is proportional to yo we would expect B = 1 and C = 3. 
The integrated SZ flux Y {= J ydVt) scales with mass as (Barbosa et al. 1996) 

Y o^M^/^{l + z)/DA{zf , (03) 

where Q is the solid angle subtended by a cluster with radius R at redshift z and Da{z) is the angular diameter 

distance to redshift z. 

For unresolved detections, where the signal to noise is simply proportional to Y divided by the noise per resolution 
element, we would expect B = 5/3. The redshift scaling is more complicated: at z ~ 0.6, oc (1 + z)^-^ leading to 
C~ -1.6. 

For resolved detections, the cluster is spread over several resolution elements: the number of detector resolution 
elements subtended by a cluster of size R scales as R^ /D\ and thus the integrated noise Ni^t scales as M^/^{1 + 
z)~^Da{z)~^ (assuming white noise). If C, traces the integrated quantity Y/Nint, we obtain ( oc M^/^(l + 2;)^d^^ and 
thus would expect i? ^ 4/3 and C ~ 0.7 at z 0.6. 

Prom these scaling arguments we would expect the slope parameter B to reside within the range 1 < -B < 1.66, 
while C is strongly dependent on the experiment. In practice, the true scaling is likely to be influenced by a number 
of observational effects. For instance, small clusters are partially washed out by the 1' SPT beam. On the other hand, 
large, low-redshift clusters are de-weighted by the matched filter due to 0MB confusion, reducing their detection 
significance (the latter could be alleviated by multi- frequency cluster detection). 

To better capture the possible form of the redshift evolution in the SPT mass-significance scaling relation, a second 
parameter could be incorporated (for example, an additional power-law dependence on the angular diameter distance 
or linear evolution in the parameter C with redshift). Given the small size of the SPT sample and the uncertainties 
placed on the other parameters, the statistical power is insufficient to constrain such a parameter. X-ray follow up 
observations of the SPT cluster sample are currently underway and will shed much light on the scaling between SPT 
detection significance and SZ flux. 

During the cosmological analysis (§5), the parameters A, B, C, and the intrinsic scatter were marginalized over 
with the Gaussian priors given in Table 01. The preferred scaling relation parameter values for the highest likelihood 
point in each of the four MOMO chains are given in Table 01. The recovered values and uncertainties on erg and w 
are unaffected by widely varying priors on both B and C, indicating that these are not a significant source of overall 
uncertainty in the current analysis. 



TABLE CI 
Preferred Scaling Relation Parameters 



Chain 


A 


B 


C 


Scatter 


ACDM WMAP7+SPT 


5.62 


1.43 


1.40 


0.21 


ACDM CMBall+SPT 


5.46 


1.43 


1.38 


0.21 


wCDM WMAP7-FSPT 


6.01 


1.42 


1.45 


0.21 


wCDM WMAP7+BAO+SNe+SPT 


5.42 


1.44 


1.41 


0.21 


(Gaussian Prior ± la) 


6.01 ±1.8 


1.31 ±0.26 


1.6 ±0.8 


0.21 ±0.04 



Note. — These are the values which maximize the likelihood for the highest likelihood point in each chain. For comparison, 
bottom row shows the Gaussian prior assumed on each parameter during the cosmological analysis. 
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D. MASS ESTIMATES 

To provide a mass estimate for the catalog presented here, one must estimate (. As mentioned in Section 4.8, ^ is a 

biased estimator of (^) and hence of C. In order to correct for this bias, one must make an estimate of the mass function. 
Although the shape of the mass function is in principle dependent on the cosmology, we find that the mass estimates 
are robust to exact choice of cosmology and hence assume a fiducial cosmology for the mass estimates presented in 
Table Dl. Specifically, we use the maximum likehhood point in the ACDM WMAP7+SPT chain: (Qb/i^ = 0.0229, 
nch'^ = 0.107, Ho = 72.7, r = 0.0877, 71^(0.002) = 0.966, ^(0.002) = 23.6 x 10-^°, A = 5.62, B = 1.43, C = 1.40, 
scatter = 0.21). The uncertainties in the cluster redshifts were assumed to be negligible for these purposes. 

Mass estimates were obtained for an assumed mass-^ scaling relation using Bayes' theorem and assuming that the 
prior probability for the mass is proportional to the Tinker et al. (2008) mass function: 

The probability distribution P(^|lnM) was found from combining the (logarithmic) intrinsic scatter in the scaling 
relation, as well as the (linear) scatter in ^ for each cluster. This yielded an unbiased estimate of the mass of each 
cluster and provided the statistical uncertainty on each mass. 
The uncertainty in mass estimates caused by uncertainty in the scaling relation was calculated using a first order 

Taylor expansion about the best fit parameters 6^, where bk = {A, B,C}. Changes in cluster mass estimates, m^, 
caused by a different assumed scaling relation are then given by Srrii ~ TJE^i^k ^ ^t)- The covariance of the mass 

estimates is given by cov(6mi, drrij) = ^ j^^ki ijj^ ' where cTki represents the covariance between parameters hk and 
hi- This parameter covariance was found by numerically calculating the Hessian of the 10-dimensional log-likelihood 
surface about the best fit parameters. In order to calculate the Hessian of the WMAP5 likelihood surface, we made 
use of CosmoMC package (Lewis et al. 2002). Assuming Gaussian errors, we inverted the Hessian matrix, yielding 
the parameter covariance. The systematic uncertainties in mass estimates reported in Table Dl are derived from the 
diagonal elements of cov(5mi, Sm,j). 

Uncertainty in the scaling relation parameters dominates the total uncertainty on the mass and the errors on the 
mass estimates are therefore strongly correlated. In Table Dl, the statistical and systematic errors are separated. The 
full covariance matrix is provided at http://pole.uchicago.edu/public/data/vanderlindelO/. 



TABLE Dl 

Mass Estimates for the SPT Cluster Catalog 



Object Name $ z M2oo(Pmean)(10" Mq/i"^) Msoo(PcWt)(10''' Mq/i"^) 
mass ± 68%stat ± 68%syst mass ± 68%stat ± 68%syst 



SPT- 


■CL 


J0509-5342 


6.61 


0.4626 


5, 


,09 


± 


1, 


.02 


± 


0, 


.69 


2, 


.98 


± 


0. 


.66 


± 


0. 


.40 


SPT- 


CL 


J0511-5154 


5.63 


0.74 


3. 


.49 


± 





.87 


± 


0. 


.43 


2, 


.15 


± 


0. 


.57 


± 


0. 


.26 


SPT- 


■CL 


J0516-5430 


9.42 


0.2952 


7. 


.84 


± 


1 


.29 


± 


1, 


.25 


4, 


.39 


± 


0. 


.84 


± 


0. 


.69 


SPT- 


-CL 


J0521-5104 


5.45 


0.72 


3, 


.39 


± 





.89 


± 


0, 


.39 


2, 


.08 


± 


0. 


.58 


± 


0. 


.23 


SPT- 


-CL 


J0528-5259 


5.45 


0.7649 


3, 


.31 


± 





.86 


± 


0, 


.38 


2, 


.04 


± 


0. 


.56 


± 


0. 


.23 


SPT- 


-CL 


J0533-5005 


5.59 


0.8811 


3, 


.19 


± 


0, 


.79 


± 


0, 


.37 


1, 


.99 


± 


0. 


.53 


± 


0. 


.22 


SPT- 


■CL 


J0539-5744 


5.12 


0.77 


3. 


.05 


± 





.84 


± 


0. 


.33 


1, 


.89 


± 


0. 


.55 


± 


0. 


.20 


SPT- 


■CL 


J0546-5345 


7.69 


1.0670 


4. 


.02 


± 





.72 


± 


0. 


.55 


2, 


.53 


± 


0. 


.52 


± 


0. 


.34 


SPT- 


-CL 


J0551-5709 


6.13 


0.4231 


4, 


.84 


± 


1, 


.06 


± 


0, 


.68 


2, 


.81 


± 


0. 


.67 


± 


0. 


.39 


SPT- 


-CL 


J0559-5249 


9.28 


0.6112 


6, 


.16 


± 


1, 


.02 


± 


0, 


.85 


3. 


.68 


± 


0. 


.71 


± 


0. 


.50 


SPT- 


-CL 


J2259-5617 


5.29 


0.1528 


4, 


.91 


± 


1, 


.45 


± 


0, 


.71 


2, 


.68 


± 


0. 


.82 


± 


0. 


.38 


SPT- 


■CL 


J2300-5331 


5.29 


0.29 


4. 


.37 


± 


1. 


.27 


± 


0. 


.60 


2, 


.48 


± 


0. 


.75 


± 


0. 


.33 


SPT- 


■CL 


J2301-5546 


5.19 


0.79 


3. 


.07 


± 





.84 


± 


0. 


.30 


1, 


.90 


± 


0. 


.55 


± 


0. 


.18 


SPT- 


-CL 


J2331-5051 


8.04 


0.5708 


5, 


.64 


± 





.99 


± 


0, 


.80 


3 


.36 


± 


0. 


.67 


± 


0. 


.46 


SPT- 


CL 


J2332-5358 


7.30 


0.32 


6. 


.21 


± 


1, 


.15 


± 


0. 


.94 


3. 


.52 


± 


0. 


.74 


± 


0. 


.52 


SPT- 


■CL 


J2337-5942 


14.94 


0.7815 


7. 


.86 


± 


1 


.17 


± 


1, 


.18 


4, 


.77 


± 


0. 


.86 


± 


0. 


.70 


SPT- 


-CL 


J2341-5119 


9.65 


0.9984 


5, 


.02 


± 





.82 


± 


0, 


.69 


3. 


.14 


± 


0. 


.60 


± 


0. 


.42 


SPT- 


-CL 


J2342-5411 


6.18 


1.09 


3, 


.21 


± 





.70 


± 


0, 


.43 


2, 


.03 


± 


0. 


.48 


± 


0. 


.27 


SPT- 


-CL 


J2355-5056 


5.89 


0.35 


4, 


.86 


± 


1 


.13 


± 


0, 


.70 


2, 


.79 


± 


0. 


.70 


± 


0. 


.39 


SPT- 


■CL 


J2359-5009 


6.35 


0.76 


3. 


.98 


± 





.84 


± 


0. 


.54 


2, 


.45 


± 


0. 


.57 


± 


0. 


.32 


SPT- 


■CL 


JOOOO-5748 


5.48 


0.75 


3. 


.36 


± 





.87 


± 


0. 


.41 


2, 


.07 


± 


0. 


.57 


± 


0. 


.25 



Note. — Note that the masses MsooiPcrit) (where the ovordensity is with respect to the critical density rather than the 
mean density) were calculated by converting from M2oo{pmean) assuming a Navarro- Frcnk- White density profile and the mass- 
concentration relation of Duffy et al. (2008). These masses may be underestimates for low-redshift clusters {z < 0.3), where 
a power-law scaling relation fails to fully capture the behavior of the CMB-confused SZ signal, and that given for SPT-CL 
2332-5358 should be considered a lower limit due to the known point source contamination. 



